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ABSTRACT 

We utilize the Sagittarius Window Eclipsing Extrasolar Planet Search 
(SWEEPS) HST/ACS dataset for a Deep Rapid Archival Flare Transient Search 
(DRAFTS) to constrain the flare rate toward the older stellar population in the 
Galactic bulge. During 7 days of monitoring 229,293 stars brighter than V=29.5, 
we find evidence for flaring activity in 105 stars between V=20 and V=28. We 
divided the sample into non-variable stars and variable stars whose light curves 
contain large-scale variability. The flare rate on variable stars is ~ 700 times that 
of non- variable stars, with a significant correlation between the amount of under- 
lying stellar variability and peak flare amplitude. The flare energy loss rates are 
generally higher than those of nearby well-studied single dMe flare stars. The 
distribution of proper motions is consistent with the flaring stars being at the dis- 
tance and age of the Galactic bulge. If they are single dwarfs, they span a range 
of ~ 1.0 — 0.25M Q . A majority of the flaring stars exhibit periodic photometric 
modulations with P <3d. If these are tidally locked magnetically active binary 
systems, their fraction in the bulge is enhanced by a factor of ~20 compared to 
the local value. These stars may be useful for placing constraints on the angular 
momentum evolution of cool close binary stars. Our results expand the type 
of stars studied for flares in the optical band, and suggest that future sensitive 
optical time-domain studies will have to contend with a larger sample of flaring 
stars than the M dwarf flare stars usually considered. 
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Introduction 



Stellar flares are the most dramatic example of variability seen in stars with outer con- 
vective envelopes during the bulk of their lifetimes. The study of stellar flares can inform 
several topics: understanding the observable manifestations of magnetic reconnection in stel- 
lar atmospheres; constraining the expected stellar flare rate in optical transient searches for 
rare phenomena; and exploring the importance of stellar flares on exoplanet environment 
and habitability. The study of solar white-light flares has revealed the importance of the UV 
and optical regions in the energetics of flares; although flare variations at other wavelength 
regions such as X-rays are larger, white light continuum emission domi nates the total ra- 



diated energy, by factors of roughly 100 with respect to X-ray energies ( lEmslie et al.l 12005 
Kretzschmarl 120 1 ll ) . 



The general consensus appears to be that stellar flares are produced by the same basic 
physical processes occurring in solar flares. This is despite the fact that flaring stars can 
be sign ificantly different f rom the Sun: F dwarfs ( IMullan fc Mathioudakisll2000l ). G and K 



( IQsten et al.l 120041 ). dMe flare sta rs (IQsten et al 



giants (lAvres et al.l Il999l . l2001t iTesta et al.l 12007 ) , tida lly locked RS CVn binary systems 



20051 ) . flares on hyperactive young Suns 



containing star- disk interactions (jFavata et al.ll2005l ). and very low mass dwarfs near the 
substellar limit ( jStelzer et al.l 120061 ) . At optical wavelengths, the selection for significant 
brightness increase in broad band filters means that M dwarfs are the primary observed 
flaring stars, due to the intrinsic red stellar spectrum compared with the blue flare spec- 
trum (flare continua are characterized r oughly by a blackbody of T~10 4 K compared with 
stellar T e ff <3500 K; lHawley et al.ll2003l ). dMe stars are known to undergo dramatic optical 
variations (peak enhancements in the U band can be several magnitudes over timescales of 
minutes), with mean flare energies in the U-band of about 10 32 erg — comparable to the 
largest solar flares. The flare rat es for the most ac tive dMe flare stars are 0.2-6 flares/ star/hr, 
10 4 -10 5 times that for the Sun (ILacy et al.lll976l ). Optical flare ti mescales on M dwar f flare 
stars can range from less than a minute to several te ns of minutes (lHawley et al.ll2003l ) with 
mean flare duration independent of flare amplitude (jKunkellll973l ). 



Due to their dependence on the presence and dynamics of interacting magnetic fields, 
stellar flares are a diagnostic of transient magnetic activity, in addition to more commonly 
used persistent indicators such as chromospheric/coronal emissions, starspots and rotation. 
In situations where persistent magnetic activity is lacking, transient magnetic activity pro- 
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vides a complementary diagnostic of the presence of large-scale magnetic fields. Magnetic 
activity i s a function of con vection zone depth and rotation rate, as quantified by the Rossby 
number ( INoyes et al.lll984j ). for stars that are not fully convective. This relates ultimately to 
the dynamo mechanism producing the large-scale magnetic fields whose stresses and plasma 
interactions cause the mechanical heating and reconnection observed as different magnetic 
activity phenomena. In single solar-like stars, magnetic activity is a function of age, with 
young stars exhibiting enhanc ed magnetic activity, and a consequ ent fall-off with increasing 
age accompanying spin-down ( Skumanichlll972l ; ISoderblomlll982l ). Likewise, the stellar flare 



rate appears to decline as a function of age. Results for low mass stars have suggested that 



activity timescal es for fully convective stars may be as large as 10 Gyr ( IDelfosse et al.lll998 



West et al.ll2008l ). While there have been several recent studies concentrating on the flaring 
rate in young stars (IWolk et al.ll2005t IStelzer et al.ll2000f ). the decline in activity at Gyr ages 
means that relatively short exposures of single stars will generally not reveal any flaring ac- 
tivity, and the flare rate in older stars has not been studied systematically. Statistical studies 
of stellar magnetic activity have found that the flare rate decreases with height above the 
Gal actic plane, as dynamical effects cause stars to move farther away from the plane as they 
age JWelsh et al.ll2007i iKowalski et~alll2009l : Irlilton et al.lboiObl ; IWalkowicz et al.lboilh . 



Stellar flares can also occur in binary systems. For tidally locked binaries, the tidal 
interaction forces the stars to rotate rapidly, in synchrony with the orbital period. Several 
classes of binary systems with at least one cool member are tidally locked (notably the RS 
CVn, BY Dra, and W UMa type binaries) and display enhanced chromospheric and coronal 
activity , and the RS CVn and BY Dra cases commonly demonstrate large individual flaring 
events (IGuinan fc Gimenea 119931 ). Because the tidal interaction causes the stars to rotate 
much more rapidly than if they were single, activity (including flaring) can be maintained in 
much older systems than for single stars. Algol binary systems, in which some mass transfer 
is taking place, can also display magnetic activity; as the spectral types of the cool secondary 
are similar to those seen in the RS CVn bin aries, it is understood that the magnetically active 
secondary star is the site of such activity (IGuinan fc Gimenea Il993l ). 



Recent results from searches for deep optical fast (<1 hr) transients have indicated the 
prevalence of stellar flaring and b olstered support for a better quantification of the flare rate. 
The three transi e nts id entified in iBecker et al.l (120041 ) turned out to be Galactic flare stars; 
Kulkarni fc Raul (120061 ) described this population as a foregro und "fog" potentially obscuring 
extragalactic fast optical transients (timescales of ^1000 s). iFresneau et al.l ( 120011 ) reported 
results from an astrographic plate survey at low galactic latitude, covering a large field of view 
(520 degrees 2 ) at B magnitudes between 1 and 14; 8% of stars sho wed flare events greater 



than 0.4 mag over 20-30 minute timescales. iRamsay fc Hakalal (120051 ) performed sensitive but 
shallow temporal coverage observations at intermediate galactic latitude, reaching depths of 
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V~22.5 on ~ minute timescales. With only 2 hours per field, they found 45 variable targets 
with 2 flare-like objects. From V > 16.5 mag to 22.5, ~ 0.2% of stellar sources were variable. 
This is especially problematic for sensitive transient searches like PanSTARR S and LSST, 



the lat ter possibly having as many as 50 stellar flares in each u band exposure (IHifton et al. 



2010af ). 



On normal stars, we expect temporal variations in the stellar light to come from surface 
activity, as well as from eclipses by stellar or substellar companions and variations due to 
seismic waves from the stellar interior. Flares are stochastically occurring, and tend to 
have fast timescales, although large long-duration flares can last for hours. On the Sun, 
white-ligh t flares have typica l durations of only a few minutes, with small enhancements 



(>10ppm; lKretzschmarll201 ll ) . A second source of temporal variations are due to the rotation 
of surface features on and off the stellar disk. Concentrations of surface magnetic field can 
produce intensity variations both lighter (such as faculae) and darker (such as spots) than 
the surrounding photospheric material: these features can be long-lived, producing regular 
photometric variations over mul tiple rotation p eriods if the size and location do not change 



appreciably during its lifetime ( jHussainl 120021 ) . Irregular periodic variations could be due 



to growth and evolution of these surface features. Coupled with differential rotation of 
the stellar surface, these effects add complexity to the periodic variations beyond a simple 
sinusoidal dependence. 

Because the flare rate is expected to be low on older stars, a systematic search for 
flares in an older stellar population needs a large number of stars, and involve a relatively 
long stare coupled with fast c adenc e to detect and resolve the flaring emission from any 



other variability. ISahu et al.l ( 120061 ) announced the discovery of 16 candidate extrasolar 
planets using nearly continuous 7 day Hubble Space Telescope/ Advanced Camera for Surveys 
(HST/ACS) monitoring of one field of 229,701 stars in two filters. The motivation for the 
initial study was to look for transiting extrasolar planets with lower masses and a larger range 
of metallicities in a part of the Galaxy that had not been explored previously. The color- 
magnitude diagram of the stars shows that the field is dominated by old stars associated 
with the Galactic bulge, as well as younger objects from the foreground disk population. 
The bulge population has an age estimated from isochrone fitting of ~ 10 Gyr. Due to the 
characteristics of the stellar target sample and observing strategy, this dataset is ideal for 
the purposes of studying the incidence of stellar flares in an older stellar population and 
comparing to the decline of magnetic activity seen in younger stellar populations, and thus 
forms the basis of the archival project described in this paper. The paper is laid out as 
follows: §2 summarizes the data reduction, §3 describes the analysis of the light curves, §4 
describes the results, §5 discusses the implications, and §6 summarizes. 
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Data Reduction 



The data use d in this paper originally came from the SWEEPS project, described in 
Sahu et all d2006h . The SWEEPS field (202"x202") was imaged continuously for 7 days 
in 2004 February. The observations consist of 254 exposures in V (F606W) and 265 in I 
(F814W), each with an exposure time of 339 seconds, and typical spacing between individ- 
ual observations ~8 minutes. The time -series phot o metry was obtained through the use of 
difference imaging and is described in ISahu et al.l (120061 ). For each star, a time series of 
relative intensities was constructed. Because the contrast of flares to the surrounding photo- 
sphere is highest at shorter wavelengths, we concentrate primarily on the V band time series 
in this paper. The magnitude and color of each star were used to estimate stellar parameters 
such as mass and temperature. The time series data are already calibrated and provide the 
starting point to examine evidence for variability and flares. 



3. Analysis 

The total number of stars in the SWEEPS sample brighter than V=29.5 is 229,293. 
After perusal of the light curves, a subset was retained where more than 70% of the light 
curve points have error bars no larger than three times the standard deviation. This reduced 
the number of stars to 222,657. A further cut was made on stars fainter than V=20 to select 
dwarf stars, reducing the number to 216,136. This serves as the starting point for our flare 
transient search, which we describe in more detail below. 



3.1. Sample Division and Detrending of Variable Stars 

The data format is a record of the relative flux of each star, where V re i,i is the relative 
flux of the ith time bin, defined as 

V t ^ = (AF/F) i =^^j (1) 

where Fi is the flux measured in time bin i, and F is the average flux computed from the 
entire time series for a given filter (here, the V606W filter). is defined similarly for the 
F814W filter. The quantity V^y is also expressed as AF/Fi in some figures in the paper. 
In order to characterize the light curves to search for transient variability, we determined a 
number of statistics. The first cut utilized the reduced chi-squared statistic of the relative 
flux summed over the time series. For these purposes xl is defined as 
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where Kei,i is defined above, <7j is the uncertainty associated with the measurement, and N\ c 
is the number of points in the light curve. Stars with a reduced chi-squared less than 1.5 
were termed non-variable, and those with xl >l-5 variable. The numbers of non-variable 
and variable stars using this definition are 214,181 and 1955, respectively. Figure [1] displays 
the distribution of these two classes as a function of stellar V magnitude. 

The variable star sample had underlying large-scale variations (large amplitude and long 
duration relative to single epochs), which needed to be removed in order to see any evidence 
of short-term flare variations. The following two approaches were applied to all light curves 
in the variable star sample, and the best fit was retained. In the first approach, periodic 
variations were identified using a Lomb normalized periodogram, and the light curve was 
folded to the dominant period. The light curve was then fit by a sum of up to 4 sine terms: 

maxl 

Vroi,i = ^2 a n sin(27mOi + p n )) (3) 

n=l 

where V re i,i is the relative flux at epoch i, 0j is the phase of epoch % when folded over the 
dominant period, p n is a phase offs et, and a n i s a nor malization. The number of sine terms 



maxl was set to 2 or 4, following iPojmanskil (120021 ) . An F-test determined whether the 



additional two terms produced a statistically significant decrease in the chi-squared statistic. 
In the second approach, the light curves were fit by a sum of sines with up to 8 terms, using 
the following prescription: 



max2 



V re l,i = ^2 a n Sin (27T71 + p n )) (4) 
n=l 

with t a normalized time (i, = ti/(t max — t mm )) where t max and t min are the last and first times 
in the time series, respectively, and max2 was 2, 4, 6, or 8, giving up to 16 free parameters. 
The value of maxl was taken to be larger than maxl to allow more flexibility with the 
irregularly variable stars. We performed fits to optimize the parameters. A decision as to 
which prescription (described by equation [3] or Hj) was better was made by comparing x 2 
and number of degrees of freedom v through an F-test. The fit was then subtracted from the 
original time series to detrend the light curve. We termed light curves fit better by the first 
approach regularly variable, and those described better by the second approach irregularly 
variable. A total of 1837 out of the 1955 variable objects had successful detrending. Of the 
1955 variable objects, 1443 were regularly variable and 512 irregularly variable; the number 
changed to 1335 and 502, respectively, for the variable objects with successful detrending. 
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The majority of the regularly variable objects discarded due to poor detrending had very 
short (<0.4 d) periods. The first two panels of Figure |2] illustrate the detrending process for 
a regularly variable star, while the first two panels of Figure |3] illustrate the process for an 
irregularly variable star. 



3.2. Flare Selection 



We tested a couple of different statistics in developing our algorithm to find flares. Ini- 
tially we used an na cut, filtering on individual points which were outliers from V re \ = at 
the level of na. With n >5 and a the standard deviation of a normally distributed dataset, 
Gaussian statistics indicates a low probability of such events happening, and this method 
can be used to identify outliers which may be flares. However, there are potentially a large 
number of cosmic ray hits which may show up as outliers using this method, and while it 
is theoretically possible to separate a true flare which is occurs in a single time bin from a 
cosmic ray by using the PSF shape of that exposure for that star, this method returns an 
unwieldy number of potential flares. We decided instead to employ a differ ent set of criteria 
based on a djacent t i me bin s. We define a statistic 4>vv based on the work of IWelch fc Stetson 
r ll993h and Istetsonl Jl996h . as follows: 



>vv 



V re \ 
a 



a 



(5) 



where V^y is the relative flux in epoch i (after detrending, if necessary; see §3.1), and ai is 
the error on the measurement in epoch i. This computation is done for each subsequent pair 
of temporal data, so that j — i + 1. Whereas IWelch fc Stetson! (119931 ) and IStetsonl (119961 ) 
developed this index to investigate variability in multi-band and multi-epoch photometry, the 
(fivv statistic here is applied to subsequent epoch pairs in the same band in our database. A 
similar technique was used in t he analysis of M dwarf flaring in the Sloan Digital Sky Survey 
Stripe 82 dataset discussed by IKowalski et al.l ( 120091 ). In order to limi t the n umber of false 
detections, we applied the false discovery rate analysis of iMiller et al.l (120011 ) to select flare 
candidate epochs by using a 4>yy threshold. 



False-discovery rate analysis (IMiller et al.ll200ll ; iHopkins et al.l 120021 ) allows us to set 
a critical threshold value of §vv to select flare epochs while ensuring a given percentage 
of false-positives. Th e FD R technique has bee n recently employed for flare rate studies in 
Kowalski et al.l ( 120091 ) and iHilton et al.l ( ]2010bl ). as well as other astrophysical applications. 
We begin by dividing the variable detrended (excluding those stars with poor detrending) 
star sample distribution of &vv pahs (464,761 epoch pairs) into a null distribution (non- 
flaring epochs produced from spurious deviations in the data) and a candidate distribution 
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(a mix of real flares and spurious non- flares). For the purposes of this calculation we exclude 
from consideration epoch pairs in which Q vv > but the values of V^y <0. Specifically, 
the null distribution is composed of the epoch-pairs consisting of a positive flux change 
in one epoch and a negative flux change in the other epoch (i.e., negative values of &vv'i 
207,381 §vv pairs). The null distributions for positive then negative flux changes was tallied 
separately from the distribution for negative then positive flux changes, to gauge similarity. 
The candidate distribution consists of the epoch pairs with a positive flux enhancement in 
both epochs (positive values of $yv an d V^y > 0; 130,117 <3Vv pairs). The sum of these two 
numbers is smaller than the total number of epoch pairs due to exclusion of eclipses, which 
have a positive $yy but are negative in both epochs. In Figure 0], we plot the absolute value 
of the nulls compared to the candidate distribution. The candidate distribution contains a 
significant tail that is not present in the nulls. Above our FDR-set threshold (see below), the 
number of spurious epochs divided by the number of total epochs (spurious epochs and real 
flares) is less than or equal to our pre-determined false- discovery rate, given by the variable 
a. The ratio of spurious epochs to total epochs using the null distribution comprised of 
negative then positive flux changes was 0.14, and for positive then negative flux changes it 
was 0.10. 



For a, we choose 10%, based on experience with similar datasets in iKowalski et al. 



(120091 ). Successful FDR analysis requires a good understanding (model) of the null. We fit 
a double Gaussian to model the absolute value of the null distribution for the variable star 
sample (one Gaussian for § vv <5.38, one for $yy >5.38) for the $vv range from to 111 
(the maximum value of &vv i n the candidate distribution). This was done to approximate 
the tail of the distribution which exhibited an apparent break near 0yy=5.38. From the 
double Gaussian, we calculate a p-value for each value of 3>vv in the candidate distribution, 
where the p-value for a &vv bin represents the probability that a null epoch-pair has that 
value of jjVy or gr eater. Then, given our pre-determined a, the FDR IDL prescription in 



Miller et al.l (120011 . ; Appendix B) sets the critical threshold. For a = 10%, the critical 
threshold is ~14.5. This limits the candidate distribution from 130,117 epoch pairs with 
4>vv >0 and V re i,ij >0 to 920 epoch pairs. We additionally apply the criterion that the flux 
divided by the standard deviation in the detrended light curve be > 2.5 in order to guard 
against light curves with marginal detrending. This limits our variable-star flare sample to 
229 epoch pairs for by-eye inspection. We apply the same $yy and standard deviation cuts 
to the non-variable star sample. The same distribution of null vs candidate epochs for the 
non- variable star sample is shown in Figure HI 

The bottom two panels of Figures [2] and [3] demonstrate the flare identification technique 
on variable stars, while Figure [5] demonstrates the process on a non- variable star. Out of 
214,181 stars in the non- variable sample, 17 epoch pairs survived using this rigorous filtering. 
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In the variable star sample of 1837 stars, 229 epoch pairs were retained. The next step was 
to visually inspect all candidate flare epoch pairs. The above steps were done using V^, 
and we examined the variations in I Te i at this stage as well. We eliminated suspect flares, 
examples of which include flares occurring near times of particularly poor detrending (such 
as at the peak of a longer-lasting rise in intensity) or on light curves with very large scatter in 
general. This further reduced the number of candidate flare epochs. At this stage it became 
apparent that systematically the first, 42nd and 43rd epochs were showing up as flares, and 
we removed these as well. This leaves us with 16 epoch pairs on 16 stars in the non-variable 
light curves, and 106 flare events (128 flare detections) on the variable light curve sample 
(from 89 stars), for a total of 122 flares from 105 stars. Of the 89 variable stars exhibiting 
flares, 63 had regular periodic variations and 26 were irregularly variable. Table [1] gives a 
tabular breakdown of the flare identification process described above. 



4. Results 

4.1. The Nature of the Flares 

Once the candidate flares were identified, we determined their properties. The flare 
selection described in the preceding section is primarily sensitive to the flare peaks, since 
that is where the signal is highest. In order to estimate flare properties such as total radiated 
energy and average flare luminosity, we considered points with V m \ > on either side of the 
adjacent pair identified by the (ftvv statistic to define the flare duration. By the definition of 
the <j)yv statistic (see §3.2), each flare has a minimum of two points in it. Figure [6] illustrates 
some example flares. 

We summed the intensities in the exposures identified as belonging to the flare and used 
the stellar luminosity to determine the flare energy in the F606W filter as follows: 

Ef&o&w — Ak K-cyAtj (6) 

i 

where L* is the filter-specific stellar luminosity (determined from th e color-magnitude dia - 



gram, assuming the star is at the mean bulge distance of 7.24 kpc; IClarkson et al.ll2008l ). 
Epm&w is the integrated flare energy in the F606W filter, and At, is the exposure time of the 
zth time bin in the flare, and V Te i,i is the relative flux in the F606W filter. We do not use pho- 
tometric parallaxes from the color and apparent magnitude of the stars because of concerns 
about binarity; see §4.2. The zero point of the F606W filte r was calcula ted as 2.87xl0 -9 erg 



cm 2 s 1 A 1 using the flux calibrated spectrum of Vega (jBohlinl 120071 ). Sequential V-band 



observations are 7.92 minutes apart (start time to start time), followed typically by a gap 
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due to spacecraft occupation; this leads to clumping in the histogram of elapsed times of 
the flare, going from 7.92 minutes to approximately 8 hours. The lower and upper values 
of the flare durations are 11 minutes and 33.9 minutes. Since we account only for the time 
when the star was observed in computing flare parameters, these derived quantities will be 
lower limits. The error on flare energy determination is estimated by adding in quadrature 
the statistical errors on each data point in the flare, and is typically ~15%, independent of 
stellar magnitude. The light curve errors are reduced for brighter stars, but the flare contrast 
is also smaller for these stars. No error in the stellar quiescent luminosity is included 
explicitly, as the spread in distance mo dulus of stars in the SWEEPS field is ~ ±1 kpc from 
the mean bulge distance of 7.24 kpc flClarkson et al.l 120081 ) . and any deviation of the star 
from the mean bulge isochrone (from which the luminosity is estimated) imparts additional 
uncertainties. 

We consid ered the distributions of energies of the flares detected on non-variable and 
variable stars. iKowalski et al.l (120091 ) demonstrated that the flare frequency distribution of 
a large number of stars at lower cadence is consistent with higher cadence observations of a 
single star. We calculate a cumulative flare frequency distribution by determining u, given 
the number of flares observed with a minimum energy E, where v = N/t, N being the 
number of flares (106 and 16 for variable and non-variable stars, respectively) and r the 
total monitoring time. The total monitoring time r is calculated as 



N T 



(7) 



where N s is the number of flaring stars considered in each category (89 and 16 for variable 
and non- variable samples, respectively), and T mon is the monitoring time per star (for the V 
band data only, this is 0.9848 days when adding up exposures in this filter over the course of 
the 6.96 day observation). The flare frequency distribution can then be plotted as v versus 
E, and the slope of the distribution de termined from a fit to the trend in log space, using 
the functional form log v = a + /3 log E (jLacy et al.lll976l ). The top panel of Figure [7] shows 
the flare frequency distribution of variable and non-variable flaring stars in the F606W filter. 
The derived intercepts, and slopes and errors are listed in Table [2j The slopes are the same to 
within the uncertainties, although it is apparent that the distribution for DRAFTS variable 
stars extends to higher energies (by a factor of 8) and exhibits a turnover near 10 33,5 erg. In 
order to facilitate comparison with other flare surveys, we convert from F606W to U band. 
We related the energy in the F606W band to that in the Johnson V band by the ratio of 
the FWHM of the two filters (850 A and 2324 A for V band and F606W, respectively), 
ignoring the fact that the central wav e length of the Johnson V filter lies blueward of F606W 
by ^400 A, then use the Lacy et al. (jl976 ) relation between Ey and E^ (E;y=1.8xEy) to 
estimate the energy of the flares in the U band (assuming that the same relation established 
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for flares on M stars holds for these flares as well). The bottom plot of Figure [7] shows 
a comparison between the flare frequency distribution of the variable flaring stars in the 
DRAFTS sample with the nearby cool binary YY Gem. The flare frequency distribution is 
flatter for YY Gem than for the DRAFTS stars. Table [2] lists the slopes derived from flares 
on the two sub-groups of variable stars and non-variable stars and YY Gem. 

For each flaring star, we calculate the average energy loss rate due to flares, by summing 
the radiated energy of all the detected flares for a given star, and dividing by the monitoring 
time. In general the errors on this derived quantity tracks the errors for the individual flare 
energies, except for the case of stars with multiple flares per star, where the error can be 
lo wer, about 10%. T his is plotted in Figure |8]for the DRAFTS stars, as well as the dMe stars 
in lLacy et al. The stellar bolometric l uminosity was det ermined by conversion using 

stellar V magnitude and the bul g e isoch rone of lSahu et al.l (120061) . The bolometri c lumin osity 
for the stars in the lLacy et al.l ( 1 19761 ) sample was taken from iReid fc Hawleyl ( 120051 ). An 
estimate of our detection limit is shown in blue for a flare consisting of two epochs at 3.8a 
(chosen because of the cutoff value of (fivv = 14.5, a/14.5 ~ 3.8). The deviation per epoch 
was then used to estimate the limiting r elative flux. The higher energy flares dominate the 
energy budget for high luminosity stars (lLacy et al.lll976l ). The trend of average energy loss 
rate versus stellar bolometric luminosity was computed for the ensemble of DRAFTS flaring 
stars, as well as the t wo su b-groups; these fits are also plotted in Figure [8] and tabulated in 
Table EJ lLacy et al.l (119761 ) presented a strong trend of the average flare energy loss rate of 
dMe stars with quiescent U band luminosity, whereas here we find a trend with quiescent 
bolometric luminosity. 



4.2. The Nature of the Flaring Stars 

Table IA1I in the Appendix lists the properties of the flaring stars. Of the 214,181 
stars in the non-variable sample, 16 stars flared. Of the 1837 variable stars with successful 
detrending, 89 flared, giving a flaring rate among non-variable stars and variable stars of 
0.007% and 4.8%, respectively. Thus, the frequency of flaring among the variable stars 
exceeds that among the non- variable stars by a factor of ~700. Variable stars also exhibit 
multiple flares per star (average of 1.25 flares per star), in contrast to the non-flaring stars, 
which only had one flare per star. We find that 71% of the flaring stars exhibited only one 
flare during the 7 day observing window. The remainder of the flaring stars showed multiple 
flares, with the largest at five flares per star. 

We are limited by statistics for the faint stars in our sample. As indicated in Figure HJ 
past a V magnitude of about 25, the number of stars returned as variable drops by roughly 
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a factor of two compared to the brighter stars, and the number of non-variable stars in 
this same magnitude range increases by a similar factor compared to brighter stars. If we 
consider only stars brighter than V=25, the flare rate among variable and non- variable stars 
is 87/1374 (6.3%) and 4/94,053 (0.004%), respectively, and the flaring frequency amongst 
these brighter variable stars exceeds that of similarly bright non-variable stars by a factor 
of 1600. The non-variable stars tend to lie at the faint end, which suggests that the signal- 
to-noise of the data are limiting our ability to see underlying variations. These stars may 
be exhibiting variability at the same relative level as the brighter stars, but our ability to 
diagnose this is compromised due to statistics. Figure M displays the distribution of the 
relative flux at the peak of the flare against V magnitude. Flares with an increase in flux 
of a few percent to a few tens of percent can be found for the stars at the bright end of the 
sample, and the flux increase required for a flare identification becomes larger as the stars 
become intrinsically fainter. This is a reflection of the statistics needed to identify flares 
and other longer temporal trends, which are better constrained for the brighter stars in the 
sample. 

There is a marked increase in the frequency of flaring amongst photometrically variable 
stars compared to the non-variable stars in our sample. Of the flaring stars brighter than 
V=25, 96% exhibit some kind of variability, as gauged by the xl statistic. The bottom panel 
of Figure [9] quantifies this further, using the "range" of the flare light curve, taken as the 
difference between the values in the 95th and 5th percentile in the original li ght curve (before 



anydetrending has taken place), and the peak flare amplitude for that star ( jWalkowicz et al. 



20111 ). For stars with multiple flares, this represents the largest amplitude flare seen. The 
non-variable stars have both the highest range and the largest flare amplitude, both likely 
reflecting the lower signal-to-noise of these generally fainter objects. 

Figure [TU] displays the locus of the flaring stars on a color-magnitude diagram, along with 
the other stars in the SWEEPS field. In general the flaring stars lie near the bulge isochrone, 
but there are several outliers, and the flaring variable stars tend to lie systematically to the 
right of the bulge isochrone. This could be due to one of three factors: (1) young disk stars; 
(2) binary stars within the bulge; or (3) stars with an enhanced metallicity. The bulge field 
analyzed in this paper lies at low galactic latitude (£>=— 2.65°), and so we may encounter a 
substantial fraction of young disk stars. Note that we cannot exclude an enhanced metallicity 
as an explanation for the offset of the flaring stars from the bulge isochrone, but we can 
examine the other two possibilities. 



We make use of the results of IClarkson et al.l (120081 ). who computed stellar proper 



motions for the stars in this field using a second epoch of data in 2006. They were able 
to measure proper motions for more than 180,000 objects, and kinematically separate bulge 
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stars from disk stars. They determined that the fraction of disk stars in the SWEEPS field 
is approximately 14% for stars brighter than I~24. Of the 105 flaring stars identified in the 
current paper, proper motions are available for 83/89 variable stars and 14/16 non- variable 
stars. Figure [TT] shows the distribution of proper motions for the flaring stars, as well as a 
random distribution of 20,000 stars in the SWEEPS field. The flare star distribution is very 
similar to the general distribution of stars in the SWEEPS field, with no preference for being 
in the disk. The proper motion constraints are consistent with the flaring stars being at the 
distance of the bulge. 

Because a large fraction of the flares detected occurred on variable stars (89 out of 105 
stars), we investigated further the properties of the underlying variability. Of the 89 variable 
stars, 63 exhibited regular periodic variations. We note that the periods are returned from 
an automatic fitting of the light curves; they correspond to the peak frequency of the Lomb 
Normalized Periodogram and are used in a utilitarian manner in §3.1 to detrend the light 
curves. This analysis doesn't examine whether the period is significant, nor does it attempt 
to address biases which can creep into periodogram analyses (e.g. from red noise or period 
aliasing). Also note that irregularly variable stars may have photometric periods which are 
statistically significant; the detrending process in section 3.1 concentrates on which of the 
two methods (given by equations 3 and 4, respectively) can better detrend the light curves. 
A plot of the period distribution of the flaring stars which exhibit regular variations shows 
that 53 out of 63 have periods < 3 days (Figure [T2l . This is ~ 1/2 the monitoring time 
of 6.96 days and may explain the apparent drop-off in the period distribution past 3 days. 
We note that three objects have apparent photometric periods in excess of the 6.96 day 
monitoring time, and are likely spurious. We concentrate on the objects having photometric 
periods less than 3 days as these are likely the most reliable, but they are subject to the 
caveats described above. Regular modulation at these timescales (for those less than about 
3 days) may be due to rapid rotation of spotted young single stars. However, the proper 
motion constraints appear to rule this out. Another possibility is that of magnetically active 
(spotted) binaries (e.g. RS CVn, BY Dra, or W UMa binaries). Most of the signposts of 
youth (fast rotation, enhanced magnetic activity as evidenced by Ha emission, X-ray emis- 
sion, starspots, flaring) are also present in active binaries, as the tidal locking enforces fast 
rotation and consequent magnetic activity signatures at old age. Note that active binaries 
are typically detected via their enhanced chromospher ic and /or X-ray emission in old stellar 



popu l ations in globular clu sters and old open clusters ( iBassa et al.ll2004l ; Ivan den Berg et al. 



2004 ; iKashyap et al.l 120061 ). so if tidally locked, the stars could still be at t he distance and 



age o f the bulge population. As the binary fraction in the bulge is uncertain ( jClarkson et al. 



20081 ). and the proper motion measurements seem to indicate that the flaring stars are bulge 



members, the scenario that these flaring stars are detached magnetically active binaries is 
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likely. Spectroscopic observations would be needed to confirm or refute the binary nature of 
the flaring stars in favor of youth. 

We can esti mate the number of a ctive binaries ex pected in this field fr om the disk, using 
the formalism of iTinney et all (119931 ). as discussed in iKoenig et all ( 120081 ). We calculate an 
effective volume V e fr which takes into account spatial inhomogeneities probed by our sight- 
line, using an exponential fall-off with increasing distance from the galactic plane. This V e ff is 

Vea = Q {^\b\) [ 2 -(£ 2 + 2 £ + 2 )* e ^(-0], ( 8 ) 

where £ = cisin \b\/h, h is the exponential scale height, \b\ the absolute value of the galactic 
latitude of the pointing (here 2.65°), Q the solid angle subtended by the 202" x 202" arcsec 
field of vi ew of the ACS obser vations, and d is the mean-bulge distance, taken to be 7.24 
kpc (from IClarkson et all 120081 ) . We use the sp ace density of actiy e binaries from near- field 
measurements, with uab =3.7x10 -5 stars pc -3 (IFavata et al.lll995l ) to determine the number 
of active binaries from the disk we would expect to see, using Nab = tt-abKa:, and find that 
Nab is between 1 and 2. This calculation shows that the expected number of disk active 
binaries is low, and cannot account for the total number of flaring objects, consistent with 
other constraints suggesting the bulge nature of the DRAFTS binaries. 

Under the assumption that the variable flaring stars in our sample are all active binaries, 
we can take a uniform distribution through the volume of the bulge sampled i n this exposure 



to con strain the space density of active binaries in the bulge. According to IClarkson et al. 



( 120081 ) the spread in distance modulus of the SWEEPS field is about 1 kpc from the mean 
bulge distance of 7.24 kpc. For a volume using the field of view translated into distance at 
the near and far end of this distance spread and height of 2 kpc, we calculate a space density 
of active binaries of 8. 8xl0~ 4 stars pc~ 3 , or an enhancement of roughly 20 compared to the 



value found locally by IFavata et all (119951 ) . The disparity between the elapsed time covered 
in the study (6.95 days) and the monitoring time per star in the F606W band (0.98 days) 
may have led us to miss flaring behavior on other variable stars if the flares occurred during 
occultations or I814W observations. The number of variable stars is roughly a factor of 20 
larger than the subset of flaring variable stars, so this explanation would not account for all 
variable stars having flares within a seven day period. It does suggest that the enhancement 
of roughly 20 in the space density of active binaries in the bulge compared to the local value 
may be a lower limit. Stellar densities are expected to be higher in the bulge compared to 
the value locally, and this may lead to an enhancement in the binary fraction. The study 
presented in this paper, while not specifically tuned to finding binaries, may be an indirect 
way to detect the presence of magnetically active binaries through sensitive flare searches 
and thus place a quantitative constraint on binary fraction in stellar associations like the 
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Galactic bulge and old open clusters. 



Discussion 



The sample of stars surveyed for flares in this paper is markedly different from the 
types of stars and wavelength ranges usually surveyed for stellar flares. Because of the large 
amount of contrast with underlying coronal emission and the impulsive nature of flares, 
many s tudies of flares haye concentrated on the high energy (extreme ultraviolet to X-r ay) 
region JWolk et all l2005[ Istelzer et all I2OO0I : kudard et allboool : bsten fc Brow^l Il999h to 
characterize flares on young stars, active stars and binaries. Due to the intr insically red light 



of M stars in quiescence, their blue flare s are easily studied in the optical ( jLacy et al .11 1976 



Kowalski et al.ll2009l ; iHilton et al.ll2010bl ). More recently, the long stellar time series returned 
by the Kepler mission h ave been explored in th e context of white light flares on stars later 
than KO spectral type (IWalkowicz et al.l 1201 ll ) for long duration flares lasting more than 
1.5 hours. The high energy measurements have been flux- and therefore distance-limited, 
and the optical flare studies have preferentially focused on younger stars, due to the known 
decline of magnetic activity and flaring with age (and the correlation between stellar age and 
increasing distance from the Galactic plane). This study is the first to cover a wide range 
of cool stellar spectral types at optical wavelengths, while sampling an older population of 
stars. 

Our supposition is that our sample of flaring bulge stars is dominated by a population 
of cool close binaries displaying magnetic activity through transient flaring activity. The 
magnetic activity of such close binaries is usually studied through their coronal and chromo- 
spheric time-averaged properties; multi- wavelength flare campaigns to study the properties 
of flares on solar neighborhood binaries (jStern et al.lll992l ; lOsten et al.ll2002l 12004) have con- 
centra ted on ultraviolet, X-ray, and radio measurements. As discussed by lHenry &: Newsom 
( 119961 ). there are very few observations of optical flaring on nearby active detached binaries 
with G and K primaries - by far the most attention has been o n optical studies of late K/M 
dwarf binaries such as YY Gem ( iDoyle fc Mathioudakislll990l ). One no table excepti o n is a 
larg e optical flare on the binary HR 1099 (K1IV+G5V) observed by both lZhang et al.l (jl990l ) 
and lHenry fc Halll ( 119911 ). with peak enha ncements noted in the V b and ( at different times) 
of 0.18 and 0.42 magnitudes, respectively. iMathioudakis et al.l (119921 ) and lHenry fc Newsom 
( 119961 ) determined the optical flare rate for the nearby act ive binary II Peg ( K2 IV +dM) 
to be near 0.2 flares hr _1 using diffe r ent da ta sets, although iHenry fc Newsom! (119961 ) noted 
variation in the flare rate. iRucinskil (jl985[ ) and references therein mention that a few flares 
have been seen in contact binaries of the W UMa type (containing magnetically active stars), 
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and they appear to follow the relations of M dwarf flare stars. Such observations do not gen- 
erally suffice to compute a flare frequency distribution with energy or other flare statistics 
to which we can compare our results, and so we instead look to the statistics on optical M 
dwarf flares for comparison. 

As noted in §4.1, the flare frequency distributions of the non- variable and variable 
DRAFTS stars contain several notable differences. The DRAFTS variable flaring stars have 
both more numerous and more energetic flares than the non-variable flaring stars, with 
maximum energy roughly a factor of 10 larger. The slopes of the distributions are similar to 
within the uncertainties, however the most energetic flares on the DRAFTS variable flaring 
stars have a much steeper slope of —2.68 above 10 33 5 erg (see Table (21) The bottom panel 
of Figure [7] compares the flare frequency distribution with ene rgy of the ensemb le of stellar 
flares from DRAFTS variable stars against that determined by lLacy et al.l (11976) for optical 
flares on the M dwarf binary YY Gem. This binary, composed of two Ml stars (IStrassmeier 
20091 ) was chosen because of its general level of enhanced magnetic activity, as evidenced 
by starspots, chromospheric and coronal emission, and a high flaring rate. It also has the 
advantage of being relatively well-observed in its flare properties. The largest stellar flare 
energies from YY Gem and the DRAFTS variable flare stars are similar, assuming that the 
same conversion applies; the study of YY Gem is clearly sensitive to lower energy events. 
The distribution of flares with energy for YY Gem has an apparent break around 10 33,5 erg, 
which is also near where the distribution of flare energies in the DRAFTS variable stars shows 
evidence of a rollover. Table [2] shows the slopes of the distributions computed for this high 
energy end. The DRAFTS flaring variable stars show a much steeper trend. The few optical 
flares seen on nearby active binaries have e xhibited flare energie s that can approach 10 6 times 
the largest solar flare energies of 10 32 ergs ( iHenry fc Halllll99ll ). alth ough the evolved nature 
of the flaring binary member(s) may play a role in these extremes. iDoyle fc Mathioudakis 
(Il990l ) noted that the binary nature of YY Gem may be responsible for the change in the 
flare frequency distribution above 10 33 - 5 erg, as the increased volumes allow for more energy 
to be stored and released during flares. Figure [8] plots the average energy loss rate due 
to flares for th e flari ng stars in our sample and the sample of M dwarf flare stars studied 
( 1976T) . YY Gem has the highest average energy loss rate of the dMe stars 



by lLacy et al 



consider e d in lLacy et al.l ( 119761 ) . It even exceeds the DRAFTS stars at the same total light. 
Ferreiral (119981 ) pointed out that such energetic flares in cool binaries require strong magnetic 
fields on at least one of the stars, and possibly also a large volume. 



Lacy et al.l (119761 ) found that the luminous M dwarfs emit more energy in flares than 



the faint M dwarfs, even when one accounts for the smaller "flare visibility" on the more 
luminous stars. We find that the energy loss rate of the non-variable stars tends to be smaller 
than variable stars, even taking into account the general faintness of the former compared 
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with latter, which biases our flare detection to larger flare events (fits to the subsets and the 
entire set of DRAFTS flaring stars are given in Table [2]). This confirms the results shown in 
Figure [7] which demonstrate more frequen t and energetic flaring for the DRAFTS variable 
flare stars. The dMe stars in lLacy et all (119761 ) show a steeper trend of flare energy loss 
rate vs. luminosity than the DRAFTS flaring stars in Figure El even though the average 
energy loss rate of the DRAFTS flaring stars can be up to 30 times higher. The energy loss 
rate of these trends at the solar value suggests roughly an order of magnitude higher flare 
energy loss rate for the dMe stars than for the DRAFTS flaring stars. It is thought that 
the relation for M dwarfs turns over at spectral type dK since even the largest solar flares 
are much weaker than what an extrapolation predicts. Possibilities for the diff erent slopes 



may a rise from the different internal structures of the dMe stars considered in lLacy et al. 



( 119761 ) compared to the higher luminosity/mass DRAFTS flaring stars, or the hypothesized 
binary nature of the DRAFTS variable stars. We note that the one dMe star in that sample 
which is a close binary is YY Gem, whose average flare energy loss rate is in the range of 
the DRAFTS flaring stars. If the bulk of the DRAFTS flaring stars are indeed binaries, this 
would shift the Lboi/L sun values for DRAFTS flaring stars in Figure [8] to the left, and bring 
the values into better agreement with YY Gem. 

Overall, only a small fraction of stars in the sample flare, but their flare rates are large 
by comparison with measures of flaring in other kinds of active stars. The 105 flaring stars 
represent 0.05% of a ll the bulge stars searched for flares, taking the 14% disk fraction of 
Clarkson et all ( 120081 ) and the number of stars with V> 20. Figure [13] shows the variation of 
flaring fraction versus stellar magnitude for the stars exhibiting underyling variability. For 
each magnitude bin, th e number of flar i ng; epo chs was divided by the number of total epochs 
in that magnitude bin. IKowalski et al.l ( 120091 ) found a maximum flare rate of M dwarfs near 
the galactic plane (at b ~ —38°) for flares with Am >0.7 mag on stars with u < 2 2 of 8 flares 
hr _1 deg" 2 . The flaring stars in the Kepler field studied by IWalkowicz et al.l ( 120111 ) were 
of spectral type K0V and cooler, and had much smaller peak amplitudes, at AF/F >0.1% 
generally, but represent an approximate flaring rate of 0.03 flares hr _1 deg -2 (this number 
may actually be higher, even for the stellar range considered in the paper, since only a subset 
of stars in the 105 degree 2 Kepler field are selected for monitoring). We compute the spatial 
flare rate for the DRAFTS flaring stars using the effective monitoring time of the dataset 
(0.9848 days) and solid angle subtended by the field of view, 0.003148 deg 2 , which for 122 
flares results i n a flare rate of 1 700 fl ares hr _1 deg -2 , two orders of magnitude greater than 
that found by IKowalski et al.l (120091 ) . The bulge is of course a much more crowded field 
than any Galactic disk field, with an intrinsically larger number of stars per square degree, 
which will contribute to this enhancement. With sensitive photometry, this study is able to 
pick out much smaller flares on a larger variety of cool stars, with peak flare amplitudes as 
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small as a few percent, which may partly explain the much higher flare rate. The flare rate 
from this study can also be expressed as 1.2 flares per flaring star per day, whi ch is about 



ten tim es higher than the X-ray flare rate of 0.13 flares per flaring s tar per day IWolk et al. 
( 2005 ) determined on young single Suns in the Orion Nebula Cluster. Osten fc Brown ( 19991 ) 
determined the extreme ultraviolet flare rate among a sample of nearby active binary systems 
to be in the range 0.1-1.5 flares per day, or 0.3 flares per day taking the entire flare sample as 
one. The DRAFTS flare rate is consistent with the upper end of this distribution, suggesting 
that the high precision of this optical study may balance the higher contrast in the extreme 
ultraviolet wavelength observations. These statistics demonstrate that the flaring stars found 
in this study exhibit flaring activity at the extremes found in magnetically active stars. This 
is also supported by the frequency distributions, comparable to one of the most energetic 
nearby flare stars, as well as by the increased average flare energy loss rate compared to most 
well-studied nearby M dwarf flare stars. 

The high flaring rate has implications for current and future time-domain studies: given 
the derived spatial flare rate of 1700 flares hr _1 deg -2 , a single 15 second frame of an exposure 
with the Large Synoptic Survey Telescope (LSST) covering 9.6 deg 2 might return as many 
as 68 flares with amplitudes greater than 1% of the quiescent luminosity. This is the same 
order as the maximum of 50 flares pe r LSST exposure fro m M dwarfs at low galactic latitudes 
with flares exceeding Am >0.1 mag (IHilton et al.ll2010af ). and suggests an underestimate of 
the contribution to the variable sky from flaring cool stars. 

The existence of flaring in these old stars opens up other questions regarding the 
longevity of binary systems at this age. The traditional assumption about a decline of 



Skumanich 


1972; 


Soderblom 


1982) 



interactions, tidal effects can cause synchronous rotation, circularizing the orbit and coupling 
the stars' rotation to the orbit. This enforces fast rotation at old age, as discussed earlier. 
The synchronization time for stars with convective envelopes is relatively sh ort, approx i- 



mately 1 MY for near-unity mass ratios and orbital periods less than 3 days ( iZahnl 119771 ). 
and a binary with an orbital period < 3 days would circularize in about 1 GY. One can ask, 
though, how long such binaries can be maintained against angular momentum loss, particu- 
larly in the present case where we may be seeing flares from 10 GY binaries. The spin-orbit 
coupling which leads to enhanced magnetic activity can also result in orbital evolution of the 
system. Whereas in a single star, magnetic torques in a stellar wind lead to angular momen- 
tum loss and subsequent spin-down with time, in a close binary the coupling of the spin and 



^^More information can be found at http://www.lsst.org. 
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orbital angular momenta paradoxically leads the angular mome ntum loss from the stella r 
wind to result in orbital shrinkage and associated faster rotation ( IGuinan fc Gimenezlll993l ). 
The exact timescale is given by the binary mass ratio, the initial period, and the nature of 
the stellar wind and its effect on magnetic brakin g. The angular momentum evolution o f 
cool close binaries has been studied quantitatively ( IGuinan fc Gimenezlll993l ; IStepiehll201ll ). 
using stellar wind/velocity variations with time found for single stars. Differences in these 
two treatments, which stem from different parameterizations of the angular momentum loss, 
can lead to differences in the timescales for evolution of a cool close binary to the contact 
stage or Roche Lobe overflow of up to an order of magnitude. 

The evolution of close binaries in the bulge could also be affected by the higher number 
densities of stars in the bulge: a tertiary can affect the angular momentum evolution of clos e 
binaries through Kozai cycles accompanied by tidal friction (IFabrycky &: Tremaind 120071 ) . 
Another possibility for these stars, if spectroscopic confirmation shows that some of these are 
single, is that they are magnetically young (by virtue of their flaring and fast rotation) but 
kinematically old. They would then be si milar to some s olar n eigh borhood chromosphericall y 
young, kinematically old stars studied by lPoveda et al.l (119961 ) and lRocha-Pinto et al.l ( 120021 ). 
and a plausible explanation for their existence at such old ages would be binary coalescence 
of a low mass, short period binary system into an apparently single, rapidly rotating (and 
flaring) star. Thus, while we do not know the exact nature of the stars in our sample, their 
activity (as deduced by the flares detected on them) and their 10 GY age (from proper 
motions consistent with the bulge) coupled with the few day periods observed from the data 
suggest that they may be able to place constraints on the angular momentum evolution of 
binary stars, if they are close binaries, or signal the formation timescales and coalescence 
rates of low mass stellar mergers. 



6. Summary and Conclusions 

We repurposed the Sagittarius Window Eclipsing Extrasolar Planet Search into a Deep 
Rapid Archival Flare Transient Search. Data from the nearly 7 day stare, coupled with 
rigorous tests to select candidate flare events, allowed us to identify 122 flares on 105 stars, 
out of a total of 216,136 stars surveyed. The characteristics of the flares are similar to those 
of optical flares on well-studied, nearby flaring stars, and suggest a high flaring rate. Stars 
with underlying variability produced flares at a higher rate than stars with no discernible 
variability, by a factor of nearly 700; this factor is 1600 when considering stars between V=20 
and V=25. We interpret these underlying variations as arising from rotational modulation 
due to starspots. The flaring stars appear to be characterized by an older stellar population 



-20 - 



which is likely binary in nature, and demonstrates the capability of tidally locked systems to 
retain magnetic activity in old age. The study expands the type of stars studied for flares in 
the optical band, and suggests that future optical time-domain studies will have to contend 
with a larger sample of potential flaring stars than the M dwarf flare stars usually considered. 
A more in-depth study of the nature of the variable flaring stars in the present sample will 
provide constraints on the nature of angular momentum evolution in old cool, magnetically 
active binary stars. 
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ciation of Universities for Research in Astronomy (AURA), inc., under NASA contract 
NAS5-26555. AFK and SLH acknowledge support from NSF grant 08-07205. We thank 
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7. Appendix 

Table IA1I lists the properties of the DRAFTS flaring stars: their positions, magnitudes 
in the F606W and I814W filters, the type of variability, the period of any underlying regular 
variability, and the range of the amplitude. 



Table 1. DRAFTS Flare Search Breakdown 



Selection Criterion 


# in Variable Sampli 


3 # in Non- Variable Sample 


Stars in the field 




229,701 


Stars with V<29.5 




229,293 


Stars in which >70% of light curve bins have errors < 3 o 




222,657 


Stars with V>20 




216,136 


xl <l-5 (non-variable) or xl >l-5 (variable) 


1955 


214,181 


Variable stars with successful detrending 


1837 




Epoch pairs with <p V v >14.5, V^y >0 ( j = i + 1) 


920 


17 


Epoch pairs with V Te \ t i/(Ti > 2.5 


229 


17 


Individual flare detections 


128 


16 


# flare events 


106 


16 


Total # of flaring events 




122 


Flaring stars 


89 


16 


Total # of flaring stars 




105 
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Table 2. Flare Frequency Distributions 



Category Band Intercept Slope 





log v = 


a + (3 log E 








DRAFTS variable stars 




F606W 


46.1 


-1.43±0.14 (N= 


= 106) 


DRAFTS non-variable stars 




F606W 


57.8 


-1.81±0.45 (N= 


= 16) 


DRAFTS variable stars 




U band 


45.8 


-1.43±0.14 (N= 


= 106) 


YY Gem 




U band 


14.5 


-0.49±0.12 (N= 


=18) 


DRAFTS variable stars, ^uji 


>10 33 - 5 er, 


g 88.1 


U band 


-2.68±0.49 (N= 


=30) 


YY Gem, E Uifl >10 33 5 erg 




U band 


55.5 


-1.69±0.69 (N= 


=6) 


logE 


tot,u/-^~mon 


= a + (3 log 


Lbo\/ Lq 






DRAFTS all stars 




U band 


28.7 


0.59±0.04 (N= 


422) 


DRAFTS variable stars 




U band 


28.7 


0.53±0.05 (N= 


406) 


DRAFTS non-variable stars 




U band 


28.2 


0.33±0.04 (N= 


46) 


Lacv et al. (1976) sample 




U band 


29.6 


1.30±0.12 (N= 


:8) 
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Table Al. Properties of DRAFTS Flaring Stars 



RA 




Dec 


V 


I) 


Variability 


Period 


Range a 


(J2000) 


(J2000) 


(F606W) 


(F814W) 


Type 


(d) 




17 co cq oo 
1 ( Oo Oo.oo 


-29 


11 


52.58 


OK AQ 

zo.4y 


99 QP 

zz. yo 


nat 




n 9n 
u.ou 


17 CO C A QQ 

1 ( Oo 04. 


-29 


10 


26.85 


on QA 
zu.y4 


1 Q PI 

iy .oo 


nat 




n n9 
u.uz 


1 7 K.8 KA 1% 
1 I Oo 04.00 


-29 


13 


25.05 


99 an 
zz. yu 


91 ^1 
Zl.Ol 


nai 




n qk 
u.uo 


17 CQ f| Q7 

i i oy u.o i 


-29 


12 


6.01 


OK 1 K 
ZO. 10 


99 9n 
zo.ou 


nai 




n op 
u.zo 


17 K.8 KQ A 1 

i ( oo oy.4o 


-29 


10 


24.00 


07 QQ 

z i .yy 


OK OO 

zo.zz 


nat 




9 Kl 
Z.Ol 


17 co no 

i t oo oo.yo 


-29 


10 


27.79 


OK cn 
ZO.OU 


99 9n 
zo.zu 


nat 




n 17 
U.o ( 


1 7 K.8 ^.9. KA 
1 ( Oo 00.04 


-29 


12 


7.09 


97 KK 
Z I .OO 


OK 78 
ZO. 1 o 


nai 




1 78 
1. 1 o 


1 7 cq n U 

i i oy u.04 


-29 


11 


0.50 


OP. QP 

zo.yo 


OA 91 
Z4.0 1 


nai 




1 <\K 
l.UO 


17 CO CO 97 

1 ( Oo Oo.o ( 


-29 


11 


15.66 


OP 1A 
Z0.O4 


99 7K 
Zo. 1 


nat 




n fi9 
U.Oo 


i 7 en a 1A 
1 ( Oy 4.04 


-29 


12 


2.30 


OK P~\ 
ZO.Ol 


99 1 A 
ZO. 14 


nai 




n 99 
u.oz 


17 cn o 17 

i ( oy o.Li 


-29 


12 


18.53 


OP. PI 


ZO. 41 


nai 




n qp 
u.yo 


17 en c in 
i i oy o. ru 


-29 


13 


7.97 


oi 8n 


9n 8*? 

ZU.oO 


nai 




n 1 1 

U. 1 1 


17 CQ O QO 

i t oy z.yz 


-29 


12 


43.19 


99 1A 
ZZ.04 


91 1 1 
Zl.ll 


nat 




n n/i 

U.U4 


17 CQ /I R1 

i / oy 4.01 


-29 


13 


27.75 


OK AP 
Z0.40 


99 OA 
ZO.Z4 


nat 




n 99 
u.oz 


17 en 1 qc 

i i oy i.oo 


-29 


12 


46.19 


97 91 
Z 1 .01 


oa an 
Z4.yu 


nai 




1 01 
1 .ZO 


i i oy o.oo 


-29 


12 


37.67 


OP 91 
ZO.Ol 


99 A 9 
ZO.40 


nai 




n fii 
u.oi 


17 en n qi 

i t oy u.oi 


-29 


13 


22.41 


91 9n 

Zl.oU 


9n ns 

ZU.Uo 


var / reg 


P AA 
0.44 


n n/i 

U.U4 


1 7 K.8 KA PQ 
1 1 Oo 04. OU 


-29 


11 


36.28 


91 An 

Z1.4U 


9n 1 7 

ZU. 1 ( 


vai / i eg 


n 1Q 
U. 10 


n n7 

U.U 1 


17 K% K7 aq 

1 ( Oo ( .4U 


-29 


13 


18.26 


9n 87 
ZU.O ( 


1 Q PO 

iy .oz 


vai / i eg 


1 AQ 

i .4y 


n 1 n 

U. 1U 


17^8 kq nn 
i i oo oy.uu 


-29 


11 


4.57 


9n QA 

zu.y4 


1 Q PI 

iy.00 


vai / i eg 


1 9^. 

1 .00 


n qp 
u.uo 


1 7 co co n7 

1 ( Oo Oo.U ( 


-29 


12 


13.00 


91 88 
Zl.oo 


9n K7 
ZU.O ( 


vai / i eg 


1 Q9 
i .yz 


n i a 

U. 14 


17 58 54.06 


-29 


12 


15.28 


21.75 


20.32 


var /reg 


0.79 


0.11 


17 58 54.88 


-29 


10 


29.63 


20.22 


19.08 


var / irreg 




0.03 


17 58 54.65 


-29 


12 


11.71 


20.18 


18.94 


var /reg 


1.47 


0.05 


17 58 55.31 


-29 


11 


28.88 


20.82 


19.70 


var / irreg 




0.04 


17 58 59.25 


-29 


10 


57.39 


20.32 


18.44 


var /irreg 




0.05 


17 58 57.10 


-29 


12 


58.73 


20.99 


19.80 


var /reg 


2.02 


0.07 


17 58 56.37 


-29 


12 


54.06 


22.41 


20.77 


var /reg 


1.33 


0.06 


17 58 59.93 


-29 


11 


55.12 


20.09 


18.94 


var /irreg 




0.04 


17 58 53.74 


-29 


10 


47.26 


20.25 


18.52 


var /reg 


2.75 


0.05 
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Table Al — Continued 



RA 




Dec 


V 


I ) 


Variability 


Period 


Range a 


(J2000) 


(J2000) 


(F606W) 


(F814W) 


Type 




(d) 




1 7 59 41 

1 1 W u W .t: X 


-29 


11 


44.32 


21 93 

Zj X . Ukj 


20 53 

ZjW .utj 
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VCll / 11 1 Cli 






1 

W . X w 
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12 
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21 57 

Zj X . W 1 
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ZjW .ZjW 
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1 


.72 


07 

W . W 1 


1 7 58 55 1 4 

J. 1 «JO WW. 11 
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12 
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?1 
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1 Q Q6 


\T^V 1 TPO" 
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n D5 

W . WW 
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J. I W(J uu.U J 


-29 


12 
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ZjW . (J W 
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W . X w 
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-L I W(J Wt: . Zj W 


-29 


12 


29.68 
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1 8 86 

x o . o w 
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04 

W . Wt: 
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-29 


12 


31.97 
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20 40 

ZjW . t: W 
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Vdl / 1 


2. 


.32 


1 1 

W . X X 


1 7 58 56 88 


-29 


11 


11.69 


99 90 

Zj Zj . ZjW 


90 Q4 

ZjW . C/t: 


^raT 1 /lTTPO* 
VCll / 11 1 "Obi 






n 05 

w. w w 


1 7 58 57 ^7 

X 1 uO W I .O 1 


-29 


10 


40.02 


99 1 7 

ZjZj. X I 


90 81 

ZjW .O X 


~\rz\v / tpo" 
Vdl / l 


6 


.44 


n 1 1 

W . X X 


1 7 58 57 76 


-29 


12 


8.28 


90 45 


1 Q 97 

X C/.Zj 1 


^rfiT /ypcr 
V dl / 1 L. ^ 


1 


,50 


n 03 

w . wo 


1 7 58 56 26 

X ( WO WW. ZjW 


-29 


13 


36.83 


91 38 

Zj X . OO 


1 9 99 


VP t* /rpcr 
Vdl / 1 eg, 





.06 


04 

W . Wt: 


1 7 58 58 78 

X 1 WO WO. I O 


-29 


10 


49.51 


22 60 

ZjZj . Uv 


20 94 

ZjW . l/t: 


VP T* /TPO* 
Vdl / 1 eg, 


2 


.27 


07 

W . W 1 


1 7 58 5Q 1 7 

XI WO UJ. 1 1 


-29 


10 


21.82 


21 64 

Zj X . W A 


20 1 9 

ZjW . X U 


VP T / TPD* 

Vdl / 1 


1 


.26 


09 

w . w u 


1 7 58 55 49 

X 1 WO WW.tiZj 


-29 


13 


9.91 


94 S4 


99 "YK 

ZjZj . 


TTflT / TPD" 
V dl / 1 C - ^ 


1 


,78 


n ^5 

w .ow 


1 7 58 57 06 

X 1 tJO W 1 .WW 


-29 


13 


35.16 


20 58 

Zj W . WO 


1 9 35 


VP! T /1TTPO* 
Vdl / 11 1 C-cl 






09 

w . w u 


1 7 58 57 21 

XI WO tJ 1 .Zil 


-29 


11 


18.19 


21 32 

Zj X . OZj 


19 94 


VP T /1TTPO* 
VCll / 11 1 Ccl 






04 

W . Wt: 


1 7 58 5Q 34 

X 1 WO U J.iJt: 


-29 


10 


35.30 


21 78 

Zj X . 1 O 


20 31 

ZjW . O X 


VP! T /TPO* 
veil / i cti 





.38 


07 

W . W 1 


1 7 58 5Q 68 

X 1 WO W<_/.WO 


-29 


13 


14.44 


22 33 

ZjZj . OO 


20 65 

ZjW .WW 


VP! T / TPO* 
veil / i 





.94 


14 


1 7 58 53 72 

X I WO ?Ju. 1 Zj 


-29 


13 


28.78 


24 02 

Zjt: . WZj 


22 22 

ZjZj . ZjZj 


VP T / TPO* 
veil / i cti 





.65 


25 

W . Zj W 


1 7 58 54 Q1 

XI WO ut:. <_/ X 


-29 


10 


30.76 


20 01 

Zj W . W X 


18 73 


VP T /1TTPO* 
veil / ii i eg, 






06 

w . w w 


1 7 58 55 71 

X 1 WO WW. 1 X 


-29 


13 


30.64 


20 90 

Zj W . U W 


1 9 66 


VP! T / TPO* 
Vdl / 1 Cli 


2 


.44 


09 

w . w u 


17 58 54 44 

X 1 W O W i . j: j: 


-29 


10 


35.83 


22.37 


20.88 


VPT / TPO* 


4. 


.81 


0.10 


17 58 55.79 


-29 


10 


57.05 


23.00 


20.64 


var/reg 


2 


.84 


0.19 


17 58 54.17 


-29 


13 


21.45 


21.23 


19.99 


var / reg 


2. 


.09 


0.05 


17 58 54.90 


-29 


12 


27.10 


20.27 


19.14 


var /irreg 






0.03 


17 59 0.19 


-29 


11 


32.85 


23.06 


20.34 


var/reg 


0. 


.83 


0.11 


17 58 55.06 


-29 


13 


12.33 


23.95 


21.98 


var/reg 


2 


,01 


0.12 


17 58 59.71 


-29 


12 


17.19 


22.55 


21.05 


var/reg 


1 


.04 


0.11 


17 58 54.44 


-29 


11 


45.74 


22.34 


20.93 


var/reg 


7. 


.60 


0.08 


17 58 59.53 


-29 


10 


29.80 


20.10 


18.75 


var /irreg 






0.09 


17 58 57.20 


-29 


11 


57.55 


23.01 


21.37 


var/reg 





.38 


0.16 
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Table Al — Continued 



RA 




Dec 


V 


I ) 


Variability 


Period 


Range a 


(J2000) 


(J2000) 


(F606W) 


(F814W) 


Type 




(d) 
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0.17 


17 59 5.87 


-29 


13 


18.64 


22.24 


20.56 
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0.10 


17 59 8.41 


-29 


10 


20.53 


22.53 


21.08 


var/reg 


1 
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Table Al — Continued 



KA 




Dec 


V 


1 ) 


Variability 


Period 


Kange 


(J2000J 


(J2000) 


(r 606 Wj 


(r 814WJ 


Type 


(d) 




17 59 8.49 


-29 


12 36.43 


22.83 


21.23 


var/reg 


2.78 


0.13 


17 59 3.37 


-29 


11 29.90 


24.22 


22.30 


var/reg 


0.71 


0.19 


17 59 4.86 


-29 


12 4.81 


21.05 


19.74 


var/reg 


1.18 


0.06 


17 59 7.02 


-29 


10 55.17 


22.19 


20.62 


var/reg 


0.61 


0.15 


17 59 5.47 


-29 


12 9.13 


20.75 


19.47 


var/reg 


0.85 


0.08 


17 59 2.63 


-29 


13 34.60 


22.52 


20.99 


var /irreg 




0.09 


17 59 3.58 


-29 


10 59.49 


21.77 


20.30 


var/reg 


2.02 


0.07 


17 59 0.82 


-29 


12 30.11 


22.41 


20.92 


var/reg 


1.87 


0.11 


17 59 1.38 


-29 


12 58.27 


22.91 


21.32 


var /irreg 




0.20 


17 59 6.17 


-29 


11 59.10 


22.72 


21.02 


var/reg 


2.84 


0.10 


17 59 5.17 


-29 


10 20.90 


21.87 


20.47 


var/reg 


0.84 


0.10 


17 59 4.23 


-29 


12 52.31 


23.32 


20.48 


var/reg 


0.66 


0.13 


17 59 7.00 


-29 


13 1.54 


25.32 


24.26 


var/reg 


11.79 


0.77 


17 59 4.73 


-29 


11 58.55 


24.51 


22.38 


var/reg 


1.54 


0.29 


17 59 4.47 


-29 


11 30.23 


22.89 


21.28 


var /irreg 




0.12 



a Range is defined as the difference in amplitude between the values in the 95th 
and 5th percentile in the original light curve; see §4.2 for details. 
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Fig. 1. — Plot of distribution of stars as a function of Vp606W magnitude, for stars deemed 
variable and non- variable, respectively. See §3.1 for discussion of sample division. 
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Fig. 2. — Example of a flare identified on a star with a variable light curve exhibiting regular 
periodic variations. Top panel shows relative fluxes versus time, and second panel from top 
shows the same data points, but folded over the best-fit period; the fit is shown in red. The 
third panel from the top shows the data after subtraction of the fit; red squares indicate time 
bins identified as a flare event. The bottom panel shows the (ftw values for the detrended 
light curve. Only points which lie above the threshold value of </yy=14.5 (dotted line in 
the bottom panel), have a value more than 2.5 times the standard deviation, and exhibit a 
positive increase are identified as flare peaks. The flare event (red squares in the third panel 
from the top) is constructed from data points adjacent to the flare peaks which have relative 
fluxes greater than zero. See §3.1 and §3.2 for more discussion. 
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Fig. 3. — Example of a flare identified on a star with irregular variations. Panels are as 
in Figure [2j The irregularly variable stars may have varying amplitudes, such that when 
phase-folded over the dominant period these amplitudes don't line up, or a changing slope to 
the rise of each light curve portion. In this example, the rise in amplitude from t = 60 — 61 is 
different than from t = 63 — 65, and the peak amplitudes are different. The fitting methods 
described in §3.1 are designed to detrend the light curve so that smaller scale temporal 
variations (flares) can be identified and studied. See §3.1 and §3.2 for more discussion. 
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Fig. 4. — False discovery analysis of the variable sample (top) and non- variable sample 
(bottom). The null distributions (null samples #1 and #2) are epoch pairs with negative 
and positive flux changes (leading to negative 4>yy values - their absolute value is plotted 
here); the two null samples differentiate the ordering of the flux changes. The green curve in 
the top panel shows the sum of two Gaussians fitted to the null sample. The vertical dotted 
line indicates the cut-off value used to ensure a false discovery rate of 10%. See §3.2 for more 
discussion. 
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Fig. 5. — Light curve from a no n- variable star, showing flare peak identification and points 
associated with a flare event. The top panel shows the relative flux variations as a function 
of time; the bottom panel shows the associated value of <pvv calculated for each time bin. 
The dotted line in the bottom panel indicates the threshold above which peak flare events 
are identified. All points on either side of this identified flare peak which have relative fluxes 
greater than zero are included in the flare event. The red squares in the top panel indicate 
points identified as a flare event. 
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Fig. 6. — Close-up of five sample flares. Panels (a), (b), and (d) show flares which occurred 
on stars exhibiting underlying variability, and panel (c) depicts a flare seen on a non- variable 
star. Light curves on variable stars have been detrended, as described in §3.1. Panel (a) 
shows two flares on the same star while panel (b) shows a complex flare event. 
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Fig. 7. — ( top ) Frequency of stellar flares as a function of flare energy, for flares occurring on 
non- variable stars (open circles) and variable stars (filled circles). The computation of quan- 
tities is described in §4.1. Green lines indicate fits to the flare frequency distributions, the 
parameters of which are listed in Table |2j (bottom) Flare frequency distribution of variable 
stars after conversion to an equivalent U-band flare energy (filled circles), and compared to 
the flare frequency distr ibution of the wel l-studied cool dwarf binary YY Gem (asterisks), 
using the observations of lLacy et al.l ( 119761 ) . The red lines show fits to the distribution above 
an en ergy of 1Q 33 - 5 erg, correspond ing to a break in the distribution of flare energies for YY 
Gem ( iDoyle fc Mathioudakislll990l ) and to a roll over in the distribution of DRAFTS variable 
stars. 
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Fig. 8. — The energy loss rate due to flares plotted against stellar bolometric luminosity, for 
each star identified as a flare star. The calculation of quantities is discussed in §4.1. Filled 
black circles represent results for variable DRAFTS flare stars, while op en circles show res ults 
for no n- variable DRAFTS flare stars. Red filled circles are results from lLacy et al.l ( 119761 ) for 
a sample of M dwarf flare stars. The blue dashed line indicates the approximate limit based 
on 3.8 a detections at two epochs, which gives the cutoff value 4>w of 14.5. The deviation 
per epoch was used to get a limiting relative flux. A fit to the trend of flare energy loss rate 
vs stellar luminosity is shown for all DR AFTS stars in green, and the red line illustrates a 
fit to the dMe stars in lLacy et al.l (119761 ). extrapolated to Lboi=L sun . The axis at the top of 
the plot lists conversion between stellar bolometric light and mass using the bulge isochrone, 
and is only appropriate for single stars. The average energy loss rate for th e flarin g stars 
identified here far exceeds all of the M dwarf flare stars studied in lLacy et al.l (119761 ) except 
for the close binary YY Gem, although the steeper trend of the dMe stars would imply a 
larger energy loss rate at L bol =L sun . 
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Fig. 9. — (top) Plot of peak flare amplitude as a function of Vf606W of the star, for flares 
identified through selective filtering of the time series. Filled circles represent flares on 
variable stars, and open circles flares on non-variable stars. The dotted line represents the 
noise level given by Poisson statistics for a star of a given magnitude. In general, the statistics 
on the brighter stars enabled identification of flares with increases of a few tens of percent, 
while for the fainter stars only larger variations could be detected, (bottom) Plot of the 
peak flare amplitude versus the range of the star's light curve, defined as the difference in 
amplitude between the 95th and 5th percentiles, for all flaring stars. Symbols are as in the 
top panel. The dashed line indicates a line of equality. For stars with multiple flares, the 
peak flare amplitude is the amplitude of the largest flare. There is strong correlation between 
the range of underlying variability and flare amplitudes in this sample of stars. 
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Fig. 10. — Color-magnitude diagram (CMD) for all stars in the SWEEPS dataset. Small 
black dots give location of individual stars where density of points is low; otherwise shaded 
contours are displayed. The green line gives a bulge isochrone appropriate for a 10 GY stellar 
population, with distance modulus (m — M) =14.3, foreground extinction E(B — V) = 0.64, 
solar metallicity, and [a/Fe]=0.3. Red symbols display the locus of variable flare stars on the 
CMD, while blue circles illustrate the location of non- variable flare stars. Regularly variable 
stars are indicated with a filled circle, while the irregularly variable ones are shown with 
a hollow triangle. The axis on the right side of the figure gives the appropriate mapping 
between luminosity and mass, using the bulge isochrone in green. 
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Fig. 11. — Distribution of proper motions in galactic latitude (/i/) and longitude (fib) for 
stars in the SWEEPS field. The black dots show the locations of 20,000 random stars in 
the SWEEPS field, while the red and blue points show the values of variable and non- 
variable flaring stars , respectively, which had proper motion measurements (as described in 



Clarkson et al. 



20081 ). There does not appear to be any tendency for the flaring stars to have 



a different distribution than the bulge stars. 
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Fig. 12. — Star's Vf606W versus period for the 63 flaring stars exhibiting regular periodic 
variations. Most (84%) of these flaring stars have periods less than 3 days. The relative 
number of variable stars drops for V>24.5 (as indicated from Figured]) reflecting our inability 
to establish these underlying variations. For brighter stars there does not appear to be any 
bias in the stellar brightness versus period. The conversion betw een stellar magnit ude and 
mass is shown on the axis to the left, using the bulge isochrone of ISahu et al.l ( 120061 ); this is 
only appropriate for single stars. 
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Fig. 13. — Variation of the fraction of epochs affected by flaring against stellar magnitude, 
for all variable stars. Percentages were computed by summing the number of epochs in a 
given magnitude bin which were identified as a flare by our algorithm, and dividing by the 
total number of epochs in that magnitude bin. For fainter stars our ability to detect flares 
becomes more difficult. 



